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^^ ' ABSTRACT 

We present the discovery of photometric variability in the DQ white dwarf SDSS 
J103655. 39+652252. 2 (SDSS J1036-f6522). Time-series photometry reveals a coherent monoperi- 
k> ' odic modulation at a period of 1115.64751(67) s with an amplitude of 0.442% ± 0.024%; no other 

'j_j I periodic modulations are observed with amplitudes > 0.13%. The period, amplitude, and phase of 

C^ . this modulation are constant within errors over 16 months. The spectrum of SDSS J1036-I-6522 shows 

magnetic splitting of carbon lines, and we use Paschen-Back formalism to develop a grid of model at- 
mospheres for mixed carbon and helium atmospheres. Our models, while reliant on several simplistic 
assumptions, nevertheless match the major spectral and photometric properties of the star with a self- 
consistent set of parameters: T^s ~ 15, 500K, logg w 9, log(C/He) — —1.0, and a mean magnetic field 
strength of 3.0 ± 0.2 MG. The temperature and abundances strongly suggest that SDSS J1036-I-6522 
is a transition object between the hot, carbon-dominated DQs and the cool, He-dominated DQs. 
The variability of SDSS J1036+6522 has characteristics similar to those of the variable hot carbon- 
atmosphere white dwarfs (DQVs), however, its temperature is significantly cooler. The pulse profile 
of SDSS J1036+6522 is nearly sinusoidal, in contrast with the significantly asymmetric pulse shapes 
of the known magnetic DQVs. If the variability in SDSS J1036-f 6522 is due to the same mechanism 
as other DQVs, then the pulse shape is not a definitive diagnostic on the absence of a strong magnetic 
field in DQVs. It remains unclear whether the root cause of the variability in SDSS J1036+6522 and 
the other hot DQVs is the same. 

Subject headings: white dwarfs — stars: oscillations (including pulsations) — stars: magnetic field — 
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1. INTRODUCTION 

White dwarfs (WDs) , the final stage of stellar evolution 
for the vast majority of stars, exhibit a menagerie of spec- 
tral characteristics. This is somewhat surprising at first, 
as the extreme surface gravities of WDs should chemi- 
cally stratify thei r atmospheres on as t ronomically short 
timescales (e.g. iDupuis et al.l 119921 : iZuckerman et al.l 
I200S) . The origins of the spectral species of WDs are 
thought to be either a result of the final stages of stel- 
lar evolution, such as a late thermal pulse exha.ustin g 
the residual hydrogen in a star (e.g. Ilben et aLlll983[ ). 
or as a result of physical processes along the WD cool- 
ing track, such as convective mix i ng and accretion (e.g. 
Koested[l976t IDupuis erallll993l : pTremblav fc Bergeronl 



20081 ). An understanding of these origins is crucial for 



many areas of astrophysical interest, such as the cool- 
ing r ates (and ther efore age-dat ing) of WD populations 
e.g. iWinget et al.| [l987: Hanse n et all [200 4: von Hiopel, 



20051). the fates of planetary systems (e .g. iJural 120031 : 



Zuckerman et al.ll2007t iFarihi et al.ll2009| ). and physical 



processes during the post main-sequence evolution of 
stars (e.g.lWeideman n"1990f: lHansenlll999l : lFontaine et al.l 
[200lt IWinget fc Kep ler 200^^ 

Using spectr al data f rom t he Sloan Digital Sky Survey 
(SDSS), Lieb ert et al.l (J200l identified eighteen WDs of 
spectral class DQ whose spectra exhibit lines of neu- 
tral and/or ionized atomic carbon, while the majority 
of DQ WDs have spectral features indicative of molec- 
ular carbon, such as Swan bands. A few years later, 
another DQ WD with features of atomic carbon was 
discovered i n the field of the you ng open star cluster 
Messier 35 ([Williams et al.l l2006[ ). Both papers sug- 
gest that these WDs had helium-dominated atmospheres 
with traces of carbon sufficient to mask the helium spec- 
tral fe atures. Shortly thereafter, iDufour et al.l ([20071 
I2008al) recognized that, instead, these "hot DQ" WDs 
(Teff « 18,000 - 24,000 K) have carbon-dominated at- 
mospheres, with carbon abundances log(C/He) > 1. 

The origin of these carbon-dominated atmospheres in 
a rela tively small t e mpera, ture range is difficult to ex- 
plain. IDufour et al.l ([2008a|) propose a spectral evolution 
scenario, expanded u pon by evolutionary calculations of 
lAlthaus et al.l ([2009 ^. In this scenario, a hot WD under- 
goes a vigorous very late thermal pulse that reduces the 
residual He content of the star to ^ 10^® — 10~^Mwd. 
Immediately following the very late thermal pulse, the 
WD would look like the very hot , carb on/oxygen- rich 
WD H1504-I-65 (e.g. [Werner et al.|[2004D . The residual 
helium would then rise to the surface of the WD, caus- 
ing it to evolve to spectral type DB (helium-atmosphere), 
until a convective zone mixes this helium veneer back into 
the carbon-rich lower atmosphere at roughly the temper- 
ature where the hot DQs appear. 

Alternatively, hot DQ WDs could be remnants of the 
so-called "super- AGB" stars, stars massive enough to ig- 
nite carbon, resulting in an oxygen-n eon WD surrounded 
by a layer of carbon and oxygen; iGarcia-Berro et aTl 
([1997D predict that such a star could have no residual 
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hydrogen and little, if any, residual helium. This lat- 
ter scenario requires that all hot DQ WDs should be 
massive (> lAf©), while the former scenario allows for 
lower-mass hot DQ stars. 

WD asteroseismology is a well-established method 
for studying the parameters and interior structures of 
WDs, provided non-radial pulsations are present and 
can be detected. Simple arguments involving the ther- 
mal timescale at t he base of the convection zone led 
iMontgomerv et aD (|2008 ) to suspect that some hot DQs 
might be unstable to non-radial pulsations, and obser- 
vations inspired by this possibility led to the discovery 
of variations in the hot DQ SDSS J142625.70-h565218.4 
(hereafter SDSS J1426-1-5652). Nearly simuhaneous with 
the publication of this discovery, Fontaine et al. (200^ 
published a prediction based on detailed nonadiabatic 
calculations that hot DQs may be unstable to nonra- 
dial pulsations if sufficie nt helium is present in the at- 
mosphere. Subsequently, iCorsico et al.l ([20091 ) published 
nonadiabatic pulsation analysis of hot DQs based on 
their hot DQ evolutionary models. Crucially, all three of 
these pulsational stability analyses present significantly 
different predictions of the location of the hot DQ insta- 
bility strip in T^s and log g space. 

Observations of variable hot DQs have uncov- 
ered additional mysteries about these stars. The 
folded pulse shape of SDSS J1426+5752 differs sig- 
nificantly from tho s e of pulsating DA and DB stars 
([Montgomerv et aTl |2008| ). and the star's spectrum 
shows Zeeman splitting indicative of a str ong mag- 
netic field (~ 1.2 MG. IDufour irall[2008bl ). Follow- 
ing the discovery o f variability in SDSS J1426-I-5752, 
I Barlow et al.l ([2008) detected variability in two addi- 
tional hot DQ stars, SDSS J220029.08-b73121.5 (here- 
after SDSS J2200-0741) and SDSS J234843.30-094245.3 
(hereafter SDSS J2348-0942). They claim both vari- 
able hot DQs have pulse shapes similar to SDSS 
J1 426-I-5752. A dditional observatio ns of these two stars 
bv IDufour eFal . (2009a) agree with iBarlow et al.l ([20081 ) 
that SDSS J2200-0741 has a pulse shape similar to SDSS 
J1426-f5752, but they find that SDSS J2348-0942 has 
a pulse shape more typical of other classes of pulsat- 
ing WDs. This, combined with suggestions of Zeeman 
sphtting in the spec trum of SDSS J2200-0741, leads 
IDufour et all ((2009a) to suggest that SDSS J2200-0741 
also has a strong magnetic field, while SDSS J2348— 0942 
d oes not. 

iDunlap et al.l ([2010t ) announced the discovery of 
the DQV SDSS J133710.19-002643.5 (hereafter SDSS 
J1337-0026). This star has low-amplitude (- 0.3%) 
photometric v ariations and has a s trongly non-sinusoidal 
pulse shape ([Dunlap et al.l l2010( ). While their low- 
resolution, high signal-to-noise spectrum shows no evi- 
dence for Zeeman splitting, a higher-resolution spectrum 
shows clear Zeeman splitting (P. Dufour, in preparation). 
Could this star have a magnetic field strong enough to 
distort the pulse shape (by whatever physical mecha- 
nism), yet be weak enough to not cause noticeable Zee- 
man splitting in low-resolution spectra? Or is there no 
relation between magnetism and the p ulse shape? 

More recently, IDufour et al.l ([2011a|) discovered vari- 
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ablility in yet another hot DQ, WD 1150+012 (SDSS 
Jl 153+0056), based on variable photon flux rates in the 
FUV spectrum obtained by the Cosmic Origins Spectro- 
graph on the Hubble Space Telescope. These observa- 
tions were not long enough to establish a pulse shape for 
this object. 

Most of these studies have been conducted under the 
presumption that the variability in DQVs is due to non- 
radial pulsations. In addition to initial models sug- 
gesting th at DQVs may b e unstable to such pulsations 
(JMontgom erv et al.l [20081 : |Fontainc ct al. 2008), claims 
have been made of the detection of additional modes 
of vari ability not harmonically related to the dominant 
mode (jDufour et al.l l2009at iGreenet al.l [20091 ). If con- 
firmed, these additional modes would be evidence in 
favor of pulsations. However, other mechanisms for 
variability have been postulated, including the musings 
of M ontgomerv et all (|2008) of accretion of carbon- rich 
material in an AM CVn-like system and the hypoth- 
esis floated by iDufour et al.l ()2008b[ ) that the oscilla- 
tions could be analogous to those observed in rapidly- 
oscillating Ap stars. Other potential causes could in- 
clude short-period binary orbits or rapid rotation of these 
white dwarfs. We discuss these possibilities more thor- 
oughly in Section ??. 

The growing number of variable hot DQ WDs should 
help to address many of the conundrums surrounding 
these objects. In particular, mapping the hot DQ in- 
stability strip can test the three independent predictions 
of pulsational properties and stability. Additional pulse 
profiles can test whether the pulse shapes of variable hot 
DQs are influenced by the presence or absence of a mag- 
netic field. Finally, additional variables will help to de- 
fine further the emerging characteristics of this enigmatic 
class of variable star. 

In this paper, we announce the discovery of variability 
in another DQ star, SDSS J103655.39+652252.fl here- 
after SD SS J1036+6522. SP SS J1036+6522 was iden- 
tified bvl Liebert et al.l (I2003D as a WD with a bsorption 
lines of atomic carbon. iSchmidt et al.l (j2003f ) identify 
three sets of Zeeman-split C I lines in the star's spec- 
trum; the splitting implies a mean field strength of 2.9 
MG, or a dipole field strength of '^ 4 MG. Our spectral 
analysis (Section ??) finds that SDSS J1036+6522 is sig- 
nificantly cooler than any of the other known DQVs. In 
Section ?? we discuss the potential implications of this 
discovery for the definition and characteristics of DQV 
WDs. 

2. TIME-SERIES PHOTOMETRY AND ANALYSIS 

Positional and photometric parameters for SDSS 
J1036+6522 are taken from the Sloan Digital Sky Survey 
(SDSS) Data Release 8 (lAbazaiian et al.ll2009D . SDSS 
J1036+6522 is located at RA= 10^36™55!39, Dec= 
+65°22'52'.'2 (J2000). The star's photometric prop- 
erties via SDSS point-spread function photometry are 
g ^ 18.514 ± 0.024, u - g = -0.217 ± 0.032, and 
g — r~ —0.325 ± 0.031. The star has a significant proper 
motion, ^ig = 26 ± 1 mas yr~^; ^s = 4 ± 6 mas yr~^ 
(jMonet et al.| [2003). 

SDSS J1036+6522 was included in the iLiebert et all 



1 a.k.a. SDSS J103655.38+652252.1 inlLiebcrt ct al.l |[2003 ) and 
SDSS J103655. 38+652252.0 in Eiscnstcin ct al. (2006 ') 



(|2003D catalog of SDSS WDs with spectral features of 
atomic carbon. At the time of our first observing run, the 
spectral parameters Tcfi and log 5 were unknown. The 
star's selection was motivated by its very blue colors, 
which suggested it might be of similar temperature to 
the other known DQVs. 

Time-series photometry of SDSS J1036+6522 was ob- 
tained over several nights in late 2008, early 2009, and 
early 2010 with the Argos high-speed photometer on 
the McDonald Observatory 2.1m Otto Struve Telescope 
(Nather_& Mukadam 2004). All observations were taken 
through a 1 mm Schott glass BG40 filter. A log of ob- 
servations is presented in Table [T] 

The photometric data wer e reduced with the m eth- 
ods and pipeline described in 'Mullall y et al.l (|2005[ ) and 
improved upon in Mullally et al. (200^. Weighted aper- 
ture photometry was performed using a variety of aper- 
ture radii; the final aperture radius was chosen based 
on the lowest noise level in the resulting discrete Fourier 
transform and was typically close to the observed full- 
width at half maximum (FWHM). The light curve was 
divided by a weighted combination of one to three nearby 
reference stars chosen to be as similar in color as possi- 
ble to the target; the precise number of reference stars 
varied from observing run to observing run; fainter com- 
parison stars were excluded if their light curves were of 
low signal-to-noise due to less-than-ideal observing con- 
ditions. Differential atmospheric extinction was removed 
by fitting a second-order polynomial to the light curve, 
and sections of the light curve were excluded from fur- 
ther analysis if clearly affected by cloud, cosmic ray hits, 
obscuration by the telescope dome due to inattentive ob- 
servers, etc. The data are corrected for the motion of 
the Earth around the Solar System barycenter using the 
method of Stumpff (1980), and all UTC leap-seconds up 
to the date of observation are accounted for. No attempt 
was made at absolute photometric calibration. 

A visual inspection of the light curves reveals no ob- 
vious signs of variability in either raw or smoothed light 
curves; an example light curve from the discovery run 
is shown in the upper panel of Figure [TJ The measured 
standard deviation in the photometry of individual data 
points is 16.6 mma (1 mma = 10^^ fractional amplitude). 

In the middle panel of Figure [T] we show the discrete 
Fourier transform (DFT) of the unsmoothed data from 
the discovery (first) observing run. A clear and signifi- 
cant signal is seen at a frequency of 897.3 /iHz with an 
amplitude of 6.2 ± 1.0 mma. The DFT of the data after 
pre-whitening by this single frequency reveals no other 
significant signal (bottom panel). 

Combining data from multiple observing runs can re- 
duce the noise in the DFT and permit more precise 
measurements of the frequency, amplitude, and phase 
of the photometric modulation, if the modulations are 
coherent between each run. Our preferred methodol- 
ogy involves bootstrapping to determine a period precise 
enough to bridge multiple run s without losing a cycle 
count fe.g. IWinget et aLlll985[ ). Unfortunately, weather 
and scheduling frustrations resulted in multiple-month 
gaps between successful observing sequences that arc too 
large to bridge with this method. 

We therefore make the assumption that the photo- 
metric modulations are coherent across our entire data 
set, and we calculate a DFT for the entire combined 
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Table 1 

Observation Log 



UT Date 



Run 

Name 



Start 
(BJDtdb) 



Exp. Time 
(s) 



Run length Seeing 



(s) 



(arcsec) 



2008 Dec 29 


A1807 


2454829.88044726 


10 


2008 Dec 31 


A1811 


2454831.79331520 


30 


2009 Jan 02 


A1816 


2454833.80946591 


30 


2009 Jan 03 


A1820 


2454834.91942706 


30 


2009 Jan 04 


A1822 


2454835.78627879 


30 


2009 Apr 21 


A1858 


2454942.64579972 


15 


2009 Apr 22 


A1861 


2454943.61438557 


15 


2009 Apr 30 


A1871 


2454951.59850661 


15 


2010 Jan 16 


A2059 


2455212.82071991 


15 


2010 Feb 13 


A2076 


2455240.69756942 


15 


2010 Feb 14 


A2080 


2455241.87175010 


15 


2010 Apr 12 


A2117 


2455298.63301727 


15 



11000 


1.4 


7560 


1.8 


20430 


2.6 


11130 


1.6 


22530 


1.3 


17400 


1.1 


22455 


1.5 


15300 


1.5 


19035 


2.3 


26595 


2.2 


13758 


1.3 


17400 


1.9 



Note. — Start time is midpoint of first exposure 
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Figure 1. Time-series photometry and analysis from the discov- 
ery run of SDSS J1036-I-6522. Top panel: Photometry from 10 s 
exposures, plotted in terms of fractional amplitude f. No significant 
modulation is obvious. Middle panel: Discrete Fourier transform 
of the data in the top panel. A significant signal with amplitude 
6.2 mma is seen at 897.3 fiHz. Bottom panel: Discrete Fourier 
transform after pre-whitening the photometric data by the 897.3 
/iRz signal. No other significant peaks are present. 



data set; this DFT is shown in Figure [5] Expand- 
ing the plot about the primary peak, we see complex 
aliasing (Figure [3]) . The highest peak, with frequency 
/ = 896.340458(54) /iHz [P = 1115.64751(67) s], is the 
only alias consistent with the bootstrapped periods from 
the individual observing runs in 2008 DecembeiQ, 2009 
April, and 2010 January/February (i.e., the frequency of 
this alias is consistent with the frequency of the high- 
est peak in each of the combined runs within calculated 
errors). The window function for this period is consis- 
tent with the observed aliasing, and pre-whitening the 
light curve by this best-fit period removes all significant 
structure from the DFT (Figure [3l middle and bottom 
panels, respectively). We therefore believe that this peak 
represents the true frequency of modulation. 

We determine the amplitude and phase of the mod- 
ulation using least-squares fitting of a sine wave to the 

^ We consider all data obtained between 2008 December 29 and 
2009 January 4 as a single observing run 



Table 2 

Fits for the 1115.64751 s mode 



Run(s) 


A 


aA 


To - To^all 


(tTo 




(mma) 


(mma) 


(s) 


(s) 


Entire Set 


4.42 


0.24 


0.0 


9.5 


All 2010 


4.70 


0.39 


-1 


15 


All 2009 Apr 


3.81 


0.39 


-2 


18 


All 2008 Dec 


4.65 


0.43 


2 


16 


A1807 


5.65 


0.94 


-1 


30 


A1811 


4.10 


1.00 


-35 


45 


A1816 


3.85 


0.90 


6 


41 


A1820 


3.99 


0.80 


58 


36 


A1822 


4.42 


0.58 


-1 


23 


A1858 


5.16 


0.61 


1 


21 


A1861 


3.11 


0.65 


-45 


37 


A1871 


3.41 


0.75 


49 


39 


A2059 


5.41 


0.95 


-23 


31 


A2076 


5.02 


0.61 


-13 


22 


A2080 


4.08 


0.78 


42 


34 


A2117 


4.13 


0.81 


18 


35 



Note. — Tq — Toall is equivalent to the more fa- 
miliar O — C quantity. 

combined data set. The resulting fit gives: 

^(i)-Aosin[27r/(i-To)], (1) 

where A{t) is the observed fractional modulation, Aq — 
4.42 ± 0.24 mma is the amplitude of the modulation, 
/ is the frequency, t is the BJDtdb time of observa- 
tion (Barycentric J ulian Date of the Ba rycentric Dy- 
namical Time, e.g. lEastman et al.l l2010( ). and To = 
BJDtdb 2454830.00459(11) d. 

We also use the best-fit frequency to determine ampli- 
tudes and phases for individual observing runs. These 
results, determined again via least-squares fitting, are 
given in Table [2] and plotted in Figure |4l There is no sig- 
nificant variation of either the amplitude or To across the 
15 months, indicating that the photometric modulations 
were coherent and of constant amplitude over that time 
span. If the photometric modulations remain stable, it 
will be possible to constrain the rate of change of the pe- 
riod (P) and, potentially, discriminate between various 
models of DQVs and their evolution. 

After pre-whitening of the light curve of SDSS 
J1036-f 6522 by the best-fitting solution above, we search 
for additional pulsational frequencies. No other obvi- 
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Figure 2. Top panel: The discrete Fourier transform of tiie combination of all time-series photometry between 2008 December and 2010 
April. Bottom panel: The discrete Fourier transform of the same combined data after pre-whitening by the single significant frequency. In 
both panels, the horizontal dashed line indicates the 99% false-alarm probability level. Only one significant peak, the P = 1115.64751 s 
modulation, is observed. 
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Figure 3. Top panel: The discrete Fourier transform (DFT) for 
the combined data from 2008 December through 2010 April. Mid- 
dle panel: the DFT of the spectral window function for a single- 
frequency sinusoidal modulation with P = 1115.64751 s. Bottom 
panel: The DFT of the combined data set after pre-whitening by 
the observed modulation at P = 1115.64751 s. In each panel, 
the left portion highlights frequencies around the observed mode 
(indicated by the dashed line), and the right portion highlights fre- 
quency space surrounding the expected location of a first harmonic 
(dashed line). Pre-whitening by a single mode is removes nearly all 
power from the DFT, and any harmonic must have an amplitude 
of < 1 mma. 



Figure 4. Plots of the amplitude (top panels) and change in 
phase (equivalent to an O — C diagram, bottom panels) of the 
P = 1115.64751 s mode as a function of time. Small filled points 
indicate individual runs; open circles indicate fits to the combined 
data set in each panel. Horizontal lines are the best fit values 
(solid) and Icr errors (dashed) for the fit to the entire data set. No 
significant variability about the best-fit value is observed across the 
15-month time span. The online version of the figure uses color for 
clarity. 

ous peaks are seen in the DFT. We calculate false-alarm 
probabilities using the method of Scargle (1982) as im- 
plemented by Kepler (1993). This method gives a 99% 
false- alarm probability amplitude of 1.54 mma; no peaks 
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Figure 5. Top panel: The light curve of SDSS J1036-I-6522 folded 
at its sole observed period of 1115.64751 s (points) and a sine wave 
of the same mean amplitude and period (solid curve, red in the 
online version of the figure). Each point is the average of ~ 60 
individual photometric measurements. Two cycles are plotted for 
clarity. Bottom panel: The folded light curve after pre-whitning by 
the single observed period. No residual signal is evident; the rms 
scatter in the residuals is 1.4 mma, consistent with the photometric 
measurement errors. The residuals show no significant structure, 
in agreement with the lack of significant harmonics in the DFT. 

appear in the DFT higher than this amplitude for fre- 
quencies > 500/iHz. In partieular, no harmonic is ob- 
served at a frequency of 1793 /iHz (right panels of Fig- 
ure |3]). 

As the light curve of SDSS J1036-f 6522 is dominated 
by a single frequency of modulation, we have computed 
a pulse shape by folding the data at the 1115.64751 s pe- 
riodicity. The individual measurements are averaged in 
100 bins based on phase, giving approximately 60 pho- 
tometric points per bin. The resulting pulse shape is 
plotted in Figure [5] along with a sine wave of the same 
mean amplitude and period. Within the scatter, the 
pulse shape looks sinusoidal. 

The bottom panel of Figure [S] shows the folded light 
curve after pre-whitening the data by the single observed 
frequency. No significant residual is seen, consistent with 
the interpretation of the pulse shape as a sine wave with 
no significant harmonics. The root-mean square of the 
pre- whitened data is 1.4 mma, which is consistent with 
predictions for the scatter in each bin assuming typical 
photometric errors and Poisson statistics. 

3. SPECTROSCOPIC OBSERVATIONS AND ANALYSIS 



The SDSS spectrum of SDSS J1036-f 6522 is of mod- 
erate signal-to-noise (w 28 per resolution element at 
4600A) . This was sufficient to identify the atomic carbon 
lines a nd t he Zecman split t ing ob served by Liebcrt et al] 
(|2003| ) and lSchmidt et alJ (|2003f ). but not to obtain pre- 
cise Toff and log 5 (see Figure 6). 

In order to better constrain the atmospheric parame- 
ters of SDSS J1036-I-6522, we obtained two sets of spec- 
troscopic observations. We obtained three exposures at 
the MMT on the night of UT 2007 Dec 15 totahng 3600 
sec (3x1200) with the Blue Channel Spectrograph. We 
used the 500 line mm~^ grating with a 1" slit width, re- 
sulting in an ^ 3.6A FWHM spectral resolution over a 
wavelength range of 3400 — 6100 A. These were reduced 
using standard IRAF routines and flux-calibrated using 
the standard stars BD-f 284211 and Feige 34. 

We obtained a second data set on UT 2010 Feb 08 
with the L ow Resolution Imaging Spectrometer (LRIS, 
lOke et all llQQS) on the Keck I telescope. We took these 
observations in long-slit mode using the blue side of the 
spectrograph with the atmospheric dispersion corrector. 
We used the 600 line mm~^ / 4000 A grism in combina- 
tion with the D560 dichroic and a slit width of 1"; the re- 
sulting data have a spectral resolution of « 4.1 A FWHM. 
We obtained 2 exposures of 900 s each. 

The Keck data were reduced using the onedspec pack- 
age in IRAF. Overscan regions were used to subtract 
amplifier bias. The internal flat field was used to re- 
move pixel-to-pixel sensitivity variations; small ampli- 
tude (w 0.5%), large-scale (^ lOOA) residuals are present 
in the final flat-field. Cosmi c rays were id e ntified using 
the L.-A. Cosmic routine of Ivan DokkumI (|200lD . The 
flat-fielded exposures were averaged; pixels flagged as 
cosmic rays were exluded from the averaging. 

We extracted the one-dimensional spectrum and ap- 
plied wavelength calibrations derived from arc lamp spec- 
tra. We also derived and applied an instrumental re- 
sponse calibration derived from the spectrophotometric 
standard star G 193-74 (spectral type DC); we make no 
attempt at absolute spectrophotometry. The resulting 
spectra, shown in Figure |6l have similar signal-to-noise 
ratios of w 230 per resolution element. 

3.1. Spectral Modeling 

One major difficulty in determining physical parame- 
ters of SDSS J1036+6522 was that tools for determin- 
ing Tcff and log 5 for strongly magnetic, carbon- helium 
mixed atmospheres simply did not exist. Therefore, one 
of us (PD) set out to develop model atmospheres appro- 
priate for analyzing this star. 

3.1.1. Line Splitting 

In the presence of a weak magnetic field, typically up 
to a few Tesla (1 T = 10^ G), atoms can be described in 
the Zeeman regime. In that regime, an external magnetic 
field can be treated as a small perturbation to the Hamil- 
tonian and the energy level properties can be obtained 
from first order perturbation theory. Assuming that the 
atom is described by L-S coupling (a decent approxima- 
tion for light atoms), J and mj are good quantum num- 
bers. In the presence of a weak magnetic field, an atomic 
level with total angular momentum J splits into 2J+1 
levels with magnetic quantum number mj = -J,. . ., -|-J. 
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Figure 6. Normalized spectroscopic data for SDSS J1036+6522. 
The small differences between the continuum level of MMT (red) 
and Keck (blue) data are due to errors in flux calibration. A three- 
point average smoothing window and a downward shift has been 
applied on the the SDSS data for clarity. The majority of the 
visible features are C I with magnetic splitting. 

In that regime, the position of all magnetically shifted 
transitions between an upper and a lower level can easily 
be computed using the selection rule (Am = 0, ±1) and 
Vui = {Eu - Ei)/h, where Ei = E^q + {eh/2mec)giBmj, 
B is the field strength, Eio is the unperturbed energy 
of level i and gi is the Lande factor of level i (for re- 
cent applications in W D stars, see iFarihi et al.l 120111 : 
iZuckerman et al.|[20Tll ). 

For stronger magnetic fields able to produce splitting 
on a J-level that is comparable with the energy sepa- 
ration between the different J-levels of a term, L and 
S decouple and J is no longer a good quantum number. 
When this happens, the perturbation theory (Zeeman ef- 
fect) can no longer be applied. In this regime, called the 
Paschen-Back (PB) effect regime, the energy levels have 
to be found by the diagonalization of the total Hamilto- 
nian H = Hq + Hb where Hq is the unperturbed Hamil- 
tonian and Hb — eh/2mc (L + 2S) • B. 

Although the Paschen-Back effect was successfully in- 
terpreted within the framework of quantum mechan- 

ics long ago (e.g. see classic textbooks such as 

ICondon fc Shortlevlll935f ). very few astronomical appli- 
cations are found in non-degenerate stars. The main 
reason is that fields found in main sequence stars are 
not large enough to push most lines out of the Zeeman 
regime (sunspot field strengths are typically between 1—3 
kG, while there are cases of a few Ap stars with field up 
to 30 kG), so only when the fine-structure levels of a 
term have abnormally small separation must the partial 
Paschen-Back regime be taken into account. For such 
weak fields (by WD standards) the effect was only evi- 
denced through the very small asymmetry that it induces 
in the line profiles (,Mathvs,,1990, ). Line profile calcula- 



tions including the PB effect have been recently incor- 
porated in stellar models only for a few sensitive lines in 
Ap stars ( Fe II multiplct 74 and Li A 6707 doublet, see 
IStift et al.| [2008: Kochukhov 2008 !) and the Sun (molecu- 
lar lines and He I A10830, sec Asensio Ramos et al.' i2005l : 
IBerdvugina et aL|[2"006: Soca s-Navarr o et al. 2005) . 

Magnetic fields strong enough to push the lines into the 
PB regime ar e commonly fo und in WD stars (see the lat- 
est review bv lJordanll2009l and references therein). Yet 
since the vast majority of known magnetic WD stars only 
show hydrogen line splitting, only hydrogen-rich mag- 
netic model atmospheres have been developed in detail 
to this date. However, the population of magnetic WDs 
with heavy element feat ures is growing rapid ly thanks 
to the work on DQ (see iSchmidt et al.l [19951) , hot DQ 
(IDufour et al.ll2009bl) and DZ WDs (i.e., L HS 2534 and 
G165-7 IReid et al.ll2001t IDufour et al.l[2006D . Accurate 
determination of the atmospheric parameters for such 
stars requires the development of a new generation of 
magnetic model atmospheres that include various field 
geometries, radiative transfer for all Stokes components, 
as well as splitting in the PB regime. Although such 
models are not available yet, it is still possible to obtain 
acceptable results with a few reasonable approximations. 

3.1.2. Models and Synthetic Spectra 

As mentioned above, detailed magnetic model atmo- 
spheres appropriate for WD stars showing heavy ele- 
ments are currently unavailable. Nevertheless, it is pos- 
sible to obtain a good first order approximation with 
some reasonable assumptions. We first assume that the 
thermodynamic structure is not affected by the magnetic 
field, a conse rvative approx imation for fields of a few MG 
according to iJordanI (|1992[ ). We also calculate our syn- 
thetic spectrum using normal radiative transfer, though 
unfortunately it is currently unknown to what extent the 
use of formal radiative transfer equations for all Stokes 
parameters would affect our results. Finally, we simply 
assume a uniform magnetic field strength over the surface 
of the stars (which seems to be indicated by our obser- 
vational data, in particular the sharp and well separated 
components of the magnetically split lines in Figure 8). 
Future work, which should address all these issues, will 
be presented in due time. Equipped with this simple 
theoretical framework, we next calculate a grid of model 
atmospheres appropriate for warm DQ stars. 

Our LTE model atmosphere cod e is similar to that pre- 
sented in IDufour etHI p005) and Dufou r et al.l (j2008a| ) 
except that we use absorption line data from the Vi- 
enna Atomic Line Database (VALD). Our model grid 
covers a range of effective temperature from Toff — 
11, 500 - 19, 000 K in steps of 500 K, log g = 8.0 to 9.5 
in steps of 0.5 dex, and from log(C/He) = 0.0 to —3.5 in 
steps of 0.5 dex. 

Detailed synthetic fiuxes are then obtained from these 
atmospheric structures by first including line splitting in 
the Zeeman regime for a magnetic field of 3 MG. This 
value is obtained from measuring the separation of the 
TT (Am = 0) and a (Am = ±1) components of the well 
isolated C I 5380A triplet, shown in Figure [6l Note that 
the separation between the components for this transi- 
tion are not different in the PB regime. To estimate the 
error, we also calculated a few supplementary model with 
magnetic field slightly different then 3 MG. From this ex- 
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Figure 7. Best fit of thie Keck spectroscopic data. Insert sliows 
SDSS coverage of the He I A5876 line (a tiiree-point average 
smootiiing windowr lias been applied for clarity). Some of the 
strongest semi-forbidden transitions that were treated in the Zee- 
man regime (see text) are marked by ticks and clearly not well 
matched. Below these spectra is a synthetic spectrum with the 
same atmospheric parameters but calculated assuming splitting in 
the Zeeman regime; mismatches with the data are obvious. 

periment, we find that varying the field by more than 0.2 
MG changes the line positions sufficiently to be detected 
by visual inspection. Hence, we adopt a field strength of 
3.0 ±0.2 MG. 

While a 3 MG magnetic field is appropriate to repro- 
duce the observed splitting in the ~ 5380A features, it 
fails to reproduce the position of the carbon feature be- 
tween ^ 4700 — 5000 A (see the synthetic spectra calcu- 
lated under the Zeeman regime assumption in Figure [7]) . 
As explained above, the most likely explanation is that 
the magnetic field is sufficiently strong to push some of 
the lines in the Paschen-Back regime where the position 
and strength of spectral features differ to that predicted 
by simple Zeeman calculations. 

Indeed, as can be observed in Figure [SJ the posi- 
tion of the energy levels produced by a 3 MG mag- 
netic field can sometimes be significantly different in the 
Zeeman and Paschen-Back regime. Under such circum- 
stances, the energy levels and line strength values due to 
the magnetic field perturbation must thus be obtained 
by the diagonalization of a set of matrices for each of 
the transitions included in our synthetic spectra. We 
thus proceeded in calculating the line strength and po- 
sitions of all the components of magnetically split C I 
lines in the PB regime (the relevant equations are de- 
scribed in detail in Section 3.4 of th e monograph by 
ILandi Degl'Innocenti fc Landolfil l200l . NIST Atomic 
Spectra Database energy levels were used as input for 
the PB calculations, though the line strengths of each 
transition are obtained by re-normalizing appropriately 
the sum of all the log gf of a given multiplet taken from 
the VALD database. 

Our final grid of synthetic spectra is calculated with 
a mean magnetic field of 3 MG based on the separa- 
tion of the 5380A lines, which are transitions between 
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Figure 8. Magnetic splittings of the 3 lower sub-levels of the ^P 
term of the C I multiplet near 4775 A. Dotted lines show the energy 
level as a function of magnetic field strength assuming Zeeman 
splitting while the full line is the result of the calculations in the 
Paschen-Back regime (see text). 

^P terms that are correctly described with the Zeeman 
regime. Since the positions of the He I lines are correctly 
reproduced in the Zeeman regime, it was unnecessary to 
perform the lengthy PB calculation for helium for this 
exploratory analysis. 

Finally, we note that our C I line list contained a few 
semi-forbidden transitions that violate the selection rules 
(AS" = and/or AL = ±1,0). The PB formalism de- 
scribed above cannot be used for these lines and a more 
general approach needs to be implemented in that case. 
Unfortunately, the relevant data necessary for such cal- 
culations are currently unavailable. We thus simply mod- 
eled those transitions assuming normal Zeeman splitting 
patterns, an approximation that does not correctly rep- 
resent the positions and strength of such lines (see, for 
example, tick marks in Figure [7]). 

3.1.3. Fitting Technique, Results and Discussion 

We attempted to fit the spectra using the usual non- 
linear least square routines commonly used in spectral 
fitting, but these led to regions of parameter space that 
made little sense. For example, a fit might converge suc- 
cessfully on a minimum x^i but visual inspection would 
reveal obvious shortcomings such as not reproducing the 
447lA He I dip, not matching the slope, or matching 
the strengths and positions of one group of lines while 
another set was poorly fit. It is likely that the models 
do not perfectly reproduce reality because of the numer- 
ous approximations, and some badly reproduced regions 
(such as the group of C I lines near 4340A that are al- 
ways too strong, the semi-forbidden transitions, etc.) are 



The DQV SDSS J1036+6522 



simply pulling the solution away from realistic values. 
Even when we fixed log 5 and T^s-, the carbon abundance 
would wander to unrealistic values that resulted in obvi- 
ous visual mismatches. We therefore resorted to visual 
inspection for selecting the best model. 

Another difficulty is the strong group of C I lines near 
4340A. Choices of parameters that fit this complex well 
result in very poor fits for the rest of the spectrum; there 
is no combination of parameters in our model grid that 
can achieve a satisfactory fit over all spectral region si- 
multaneously. This could imply that our numerous ap- 
proximations affect the temperature structure in such a 
way to poorly reproduce this complex or perhaps that 
the log gf values for these lines are incorrect. 

Our visual inspection therefore starts with the C II 
lines near 4267A. Since this is an ionized feature, it pro- 
vides decent constraints on Teg. We then use the broad 
He I 4471A dip as a main indicator of carbon abun- 
dances, as this dip proves quite sensitive to log(C/He). 
In parallel, the MMT and SDSS spectra include the He I 
5876A line, which we used as a confirmation of the abun- 
dance. This line also allows some of the coolest solutions 
to be rejected. 

Our best solution, shown in Figure [71 has Toff = 15500 
K, logg = 9, and log(C/He) = -1.0. The solutions 
are degenerate; similar quality solutions can be found 
by compensating a downward change in Toff with down- 
ward changes in logg and log(C/He). We estimate our 
correlated uncertainties in these parameters as follows. 
A ^ lOOOK change in our best solution's Toff resulted in 
a ^ 0.5dex change in log(C/He), with a lower Toff re- 
sulting in less carbon and vice- versa. A change in logy 
has little effect on the carbon abundance as long as we 
stay close to logg = 9.0; at 15500K and logg = 8.5, 
log(C/He) is only 0.2 dex higher. At logg = 8.0, the 
carbon abundance increases by about 0.7 dex as com- 
pared to logg = 9. Even though logg = 9 provides the 
most appealing solution visually, given all the approxi- 
mations and the preliminary nature of our calculations, 
our value of logg should not be considered a precision 
measurement. Much more work is needed before we can 
confidently determine the mass of this object. 

Carbon-rich solutions, favored for hot DQs, do not fit 
well. However, we can impose some useful additional 
constraints on these parameters. For instance, we com- 
pared the SDSS photometry of SDSS J1036-f 6522 with 
synthetic photometry calculated by folding the model 
spectra through filter response curves. The photomet- 
ric data (Figure [9]) seem to favor a solution in the 
15000—16000 K range, a temperature range which is cor- 
roborated by the He I 5876A line, although the unknown 
reddening makes this a bit uncertain. This conclusion 
also supports the qualitative best matching model above. 

3.2. What IS SDSS J1036+6522? 

Based on the spectroscopic and photometric analysis, 
SDSS J1036-t-6522 seems to be an intermediate object 
bet ween the class i cal, c ool DQ stars (Tcff < 10, 000 K, 
e.g. IDufour et all [200I and hot DQ stars (18, 000 < 
Tcff < 24,000 K. iDufour etldl [2008aL [200l . Its high 
carbon abundance also seems to make it akin to the 
second sequence of carbon-rich cool DQs ([Dufour et al.l 
[2005; .Kocster fc Knist,,2006.) and two DQs with similar 
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Figure 9. SDSS photometric observations in the ugriz bands 
(error bars) and average model fluxes for our favored atmospheric 
parameters (circles, red in the online version of the flgure, see text). 
The photometric observations have been corrected fo r extin ction as 
prescribed by the reddening maps of lSchlegei et aQ II1998I V 



temperatu re, G35-26 (ITheill et al.l [T990l) and G227-5 
(Wcgner fc Koesteilll985D . In addition, the high mag- 
netic field of the star also appears to be a typical char- 
acteristic of many of the hot DQs (jDufour et al. 2009b) . 

Although the high apparent surface gravity should be 
considered preliminary, a growing body of evidence sug- 
gests that hot DQs and the carbon-rich sequence of cool 
DQs are massive WDs. This evidence includes unpub- 
lished parallaxes of G35— 26 and G227— 5 (Liebert et al] 
120031 C.C. Dahn, private communication); the high 
mass of th e carbon- and oxy gen-rich post-AGB star 
H1504-I-65 (jWerner et al.ll2004D : the presence of a hot 
DQ in the open star cluster Messier 3 5, in which the 
DA c l uster member WDs are all massive ([Williams et al.l 
I2OO6I [2009); and improved spectroscopic observations 
and model atmospheres of hot DQs, which favor high 
surface gravities (Dufour et al., in prep). 

We note however that these calculations are done with 
classical approximations for the Stark broadening of C I 
lines, so the high gravity solution should be considered 
highly uncertain. State-of-the-art Stark broadening pa- 
rameters for C I siin ilar to that calculated for C II in 
[Dufour et al.[ ([2011b( ) are unfortunately not available yet. 

In summary, SDSS J1036-t-6522 is most probably a 
cooler ve rsion of S DSS J1402+3818 and SDSS J0005- 
1002 (see IDufour et al. 2009b) and/or a hotter version of 
G227-5 and G35-26. SDSS J1036-I-6522 could thus be 
a transition object between the coolest carbon-rich hot 
DQs and the hottest, cool helium-rich DQ white dwarfs. 
There is thus growing evidence that this class of WDs in- 
deed forms a distinct and fairly continuous evolutionary 
sequence of mas sive hydrogen- and helium-deficient stars, 
as proposed by IDufour et al.l ([2008a) . Further analysis 
on this type of object will require improved models in- 
cluding various magnetic field geometries, solving the ra- 
diative transfer equation for all Stokes parameters (spec- 
tropolarimetric observations would be especially useful 
to test the eventual predictions from such models), and 
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inclusion of state-of-the-art Stark broadening parameters 
for C I once available. 

4. DISCUSSION 

The detection of variability in SDSS J1036+6522 in- 
creases yet again the number of known DQVs. However, 
the spectroscopic analysis raises the question of whether 
or not SDSS J1036-I-6522 belongs in the same class as 
the other known DQVs. 

There are some crucial differences between SDSS 
J1036-f 6522 and the hot DQVs, first and foremost being 
Toff, with SDSS J1036+6522 being significantly cooler 
than the other DQVs, as discussed above. Second, the 
atmospheric composition of SDSS J1036+6522 is differ- 
ent than both the hot and cool DQs. Hot DQs are 
carbon-dominate d, with log(C/He) > 0, and He usu- 
ally undetected (jDufour et al.l l200"8al) . In the "second- 
sequence" cool DQs, log( C/He) < —3 (jKoester &: Knisti 
l20M IDufour et al.l[2005h . Again, SDSS J1036+6522 is 
intermediate, with log(C/He) ~ —1. This observation 
furthers the hypothesis that SDSS J1036-I-6522 is a tran- 
sition object between hot and cool DQs. 

We now explore the interpretations and implications 
of this conclusion. 

4.1. What is the cause of the variability? 

From an observational standpoint, the photometric 
variability of SDSS J1036-I-6522 is very similar to the 
other known DQVs. These similar characteristics include 
the period and amplitude of the photometric variabil- 
ity, the dominance of a single frequency in the discrete 
Fourier transform, and the stability of the frequency and 
phase. 

Based on the time-series photometry, SDSS 
J1036-I-6522 looks and acts like the other members 
of the DQV class of objects. Naively combining this 
information with the spectroscopic evidence that this 
star belongs to the same evolutionary sequence and 
object class as the hot DQVs, we would come to the 
conclusion that SDSS J1036-I-6522 is a typical DQV and 
is variable for the same reason as the others, whatever 
that reason may be. 

The difficulties arise when we attempt to explain the 
origin of the photometric variations. As described in the 
introduction, the preferred explanation for the variabil- 
ity in DQVs has been nonradial pulsations. Even though 
there are differing physical mechanisms and different pre- 
dicted instability strips for pulsations in DQVs, SDSS 
J1036-I-6522 is far cooler than any of these predictions 
(though we note that these models do not include a sig- 
nificant magnetic field). Moreover, if our atmospheric 
analysis is correct, the atmosphere of SDSS J1036-F6522 
is helium dominated and ^ 10"^ K coo ler than the coolest 
o bserved DBV (e.g. iNjtta et al.l l2009l). 

iMontgomerv et al.l ()2008D hvpothesize that the peri- 
odic photometric variability could be due to an accret- 
ing system akin to high-accretion states of AM CVn 
systems. This speculation was based solely on the ob- 
served pulse shape of SDSS J1426-I-5752. This hy- 
pothesis now appears highly unlikely. High signal-to- 
noise spectra of DQVs show no signs of accretion (e.g. 
IDufour et al.l l2009bl l2011aD . radial- velocity variations 
have not been detected with cataclysmic- variable-like pe- 
riods (|Grcen et al.,,2009^ although orbital periods com- 



mensurate with the photometric period have not yet been 
ruled out), and no low- accretion state systems analogues 
have been identified. Given that SDSS J1036-K6522 has 
a sinusoidal pulse shape, we find no evidence that SDSS 
J1036-I-6522 is an interacting binary systern^ 

Anot her possibility, discussed by IDufour et al.l 
(l2008b[ ). is that this star is a WD analogu e of rapidly 
oscillating Ap (roAp) stars ('e.g. lKurt3ll982D . Given the 
presence of a strong magnetic field, this oblique pulsator 
model seems to be a distinct possibility. 

It is also plausible that SDSS J1036+6522 is a rapidly 
rotating WD, with the photometric variations due to the 
presence of a persistent spot, although several arguments 
could be made against this interpretation. Accreting 
magnetic objects such as AM Her systems show photo- 
metric modulation due in part to the rotation of the WD 
carrying the accretion column in and out of view, but 
SDSS J1036-I-6522 shows no evidence of accretion. In 
addition, the measured rotation periods of single mag- 
netic WDs are on the order of hou rs and days to cen- 
turies (jSchmidt fc Norsworthvlfl99l|) , significantly longer 
than the photometric period of SDSS J1036-F6522. As 
the measured rotation rates of other classes of isolated 
puls ating WDs are also on the order of hours to days 
(e.g. iWinget fc Keple^l2008t iFontaine fc Brassardll2008l 
and references therein), the rotation hypothesis would 
require SDSS J1036-I-6522 to be an outlier among single 
WDs. Time-resolved spectropolarimetry may be able to 
determine if SDSS J1036-1-6522 is rapidly rotating. 

At present, it is unclear what the cause of the photo- 
metric variability in SDSS J1036-(-6522 is, or even if the 
variability has the same cause as for the hot DQVs. 

4.2. Magnetic fields and pulse shapes 

Typically, large amplitude DAV and DBV pulsators 
tend to have relatively non-linear (non-sinusoidal) pulse 
shapes, while low-amplitude pulsators tend to be more 
sinusoidal. This relation does not appear to hold 
among the known DQVs with sufficient data to cal- 
culate pulse shapes. In the prototype DQV SDSS 
J1426-I-5752, the amplitude of the first harmonic (6.7 
mma) is less than half that of the fundamental (17.5 
mma; IMontgomerv et al]|2008f) , while in the much lower 
amplitude DQV SDSS J2200-0741, the first harmonic 
is similar in amplitude to the fundamental frequenc y 
at w 7 mma (jBarlow et al.l 120081 : IDufour et aDl2009al l. 
These significant harmonics and their phases result in 
decidedly non-sinusoidal, boxy pulse shapes with long 
brightness maxima and short, narrow minima, the op- 
posite to the pulse sh apes found in the high ampli- 
tude DAVs and DBVs (jMon tgomerv et al."2008'). How- 
ever, t he similar-amplitud e DQV SDSS J2348-0943 (w 8 
mma; iBarlow et alTl2008t JDufour et al.l l2009all is sinu- 
soidal, an d the low-amplitud e DQV SDSS 1337-0026 
(« 4 mma lDunlap et al.ll2010D has a very non-sinusoidal 
pulse shape. 

Further, the three DQVs with non-sinusoidal pulse 
shapes all have detectable magnetic fields, while the si- 
nusoidal SDSS J 2348-0943 d oes not. This has led to 
the suggestion bv lGreen et al.l (J2009) that the magnetic 
field may be responsible for the non-linearities observed 
in many DQVs. 

The folded pulse shape of SDSS J1036-h6522 appears 
to be sinusoidal. This qualitative statement is reinforced 
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quantitatively by the lack of any observed harmonics in 
the DFT. As SDSS J1036+6522 is strongly magnetic yet 
has a sinusoidal pulse shape (within the limits of exist- 
ing data), it is clear that strong magnetic fields in DQVs 
do not always cause boxy pulse shapes. Strong magnetic 
fields in SDSS J1337-0026, SDSS J1426+5752 and SDSS 
J2200— 0741 may well be the cause of their unique pulse 
shapes, but SDSS J1036-I-6522 shows that a sinusoidal 
pulse shape in a DQV does not indicate the absence 
of a strong magnetic field, at least assuming that the 
mechanism behind the variability is the same for SDSS 
J1036+6522 and the hotter DQVs. 

4.3. Non-variable hot DQs 

iMontgomerv et all ()2008[ ) claim that variability was 
not detected in several hot DQs, as well as that those 
stars would not be expected to vary given their Tcs 
and log 5 and given predictions for an instability strip. 
In light of the low amplitude of variability in SDSS 
J1036+6522 and SDSS J1337-0026 (|Dunlap et al.l[2010l) 
and the typical detectio n threshold of ^ 5 — 10 mma 
in the data obtained for IMontgomerv et all (2008), we 
feel compelled to mention that many of the stars not ob- 
served to vary could in fact be variable at the same level 
as SDSS J1036-I-6522. Photometric precision of equal to 
or better than 1 mma would be required to rule out SDSS 
J1036-H6522-like variability. 

Therefore, we feel that the claims o f non-detection of 
variability in IMontgomerv et al.l (J2008D should be viewed 
wi th the caveat that the se limits are not very significant. 
As lDufour et al.l (|2011aD detect UV variability in four out 
of five target hot DQs, it could be that most or all hot 
DQs are in fact DQVs. A re-analysis of the time-series 
data for stars not observed to vary is underway, and more 
stringent limits on non-detections of variability will be 
presented in a forthcoming paper. These limits will also 
be crucial for testing competing models for these stars, 
such as the different proposed instability strips discussed 
in the introduction. 

5. CONCLUSIONS 

In summary, we present spectroscopy and time-series 
photometry of the white dwarf SDSS J1036-H6522. Our 
analyses lead us to the following conclusions: 

• SDSS J1036-I-6522 is photometrically variable with 
a period of 1115.64751(67) s and an amplitude 
of 0.442% ± 0.024%. The period, amplitude, and 
phase of the modulations are constant within mea- 
surement errors over at least a 16 month interval. 

• SDSS J1036-I-6522 has an effective temperature 
and carbon abundance located between the hot and 
cool DQs, suggesting that it is a transition object 
between those two populations of white dwarfs. 

• The star is magnetic, with a field strength of 3.0 ± 
0.2 MG. 

• The lack of observed harmonics suggest that the 
photometric modulation is nearly sinusoidal. If the 
cause of the photometric modulation is the same 
mechanism as in the hot DQVs, then the absence 
or strength of the star's magnetic field is not nec- 
essarily correlated with the pulse shape. 



• The mechanism of the photometric modulation in 
SDSS J1036-I-6522 is unclear. There is no evidence 
for accretion, the effective temperature is signif- 
icantly cooler than pulsation models predict, and 
the photometric period is significantly shorter than 
the typical isolated white dwarf rotational period. 

• Photometric variability is common among hot DQ 
white dwarfs, but the low optical amplitudes of 
many known DQVs imply that non-detections of 
variability in any hot DQ are only significant if 
photometric precision is high (^ 0.1%). 
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